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<N ■ Abstract 

■ We analyse the latest HI large rapidity gap data to obtain diffractive parton distribu- 

tions, using a procedure based on perturbative QCD, and compare them with distributions 
[- — , obtained from the simplified Regge factorisation type of analysis. The diffractive parton 

densities and structure functions are made publically available. 

G\ 

The HI Collaboration have recently released measurements of diffractive deep-inelastic scat- 
tering (DDIS), 7*p — > X + p, with more precision than hitherto. They have presented data 
^ ; collected using their forward proton spectrometer (FPS) [1] and also data selected using the 
large rapidity gap (LRG) method [2]. In addition they have used their LRG data, which has 
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much more statistics, to determine diffractive parton densities assuming that Regge factorisa- 
tions holds. We call this a "Regge" fit@. In principle, this Regge factorisation is not justified 
theoretically, and the purpose of this note is to present diffractive parton distributions ob- 
tained from exactly the same HI LRG data [2] using a more theoretically sound procedure. It 
is convenient to call this the perturbative QCD ("pQCD") procedure, although it also contains 
non-perturbative Regge-like contributions. 

We present the pQCD procedure first. It is shown schematically in Fig. [TJ There are three 
components. First in Fig. [T](a) we have the perturbative resolved Pomeron contribution in which 
the virtualities of the t-channel partons are strongly ordered as required by DGLAP evolution: 
/ig *C • • • <C jJ 2 <C • • • <C Q 2 ■ The scale /i 2 at which the Pomeron-to-parton splitting occurs can 
vary between /i 2 , ~ 1 GeV 2 and the factorisation scale Q 2 . For fj 2 < /i 2 ,, the representation of 
the Pomeron as a perturbative parton ladder is no longer valid and instead, in the lack of a 
precise theory of non-perturbative QCD, we appeal to Regge theory where the soft Pomeron is 
a Regge pole with intercept «p(0) ~ 1.08 [5]; see Fig. (TJ^b). In addition to the resolved Pomeron 
contributions of Figs. [H(a,b), we must also account for the direct interaction of the perturbative 
Pomeron in the hard subprocess, Fig. [H(c), where there is no DGLAP evolution in the upper 
part of the diagram. 



1 The HI Collaboration call this "proton vertex factorisation" [2]. 

2 This is the conventional approach based on collinear factorisation for the scale dependence [3] and Regge 
factorisation for the xp dependence [4]. 
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Figure 1: (a) The perturbative resolved Pomeron contribution, which is the basis of the pQCD 
approach, (b) the non-perturbative resolved Pomeron diagram, which accounts for the contri- 
bution from low scales, fi 2 < /x 2 ,, and (c) the perturbative direct Pomeron contribution. 



Each diffractive DIS event is described by four variables, usually taken as x-p, f3, Q 2 and t. 
Here, xp = 1 — p' + /p + is the fraction of the proton's light-cone momentum transferred through 
the rapidity gap, (3 = x Bj /xp is the fraction of the Pomeron's light-cone momentum carried 
by the struck quark (at LO), and x Bj = Q 2 /{2p ■ q) is the Bjorken-x variable. The data are 
integrated over the squared momentum transfer t = (p — p') 2 . The variables are shown in 
Fig. [T](a). In the pQCD approach, the diffractive structure function F 2 D( ^ is given by [6,7] 

F 2 D(3) (x P ,AQ 2 )= $^2,a®a D + C *?> C 1 ) 

a=q,g ¥=G,S,GS 

where the first resolved Pomeron term corresponds to Figs. [T](a,b) while the second direct 
Pomeron term corresponds to Fig. d](c). The C^a are the same coefficient functions as in 
inclusive DIS. The diffractive^ parton distribution functions (DPDFs), aP = zq D or zg D , are 
obtained at a scale Q 2 by evolving up from /Xq with the inhomogeneous evolution equations [6,7], 

da»(xp,z,v 2 ) = ^ Paa/ fl/D + J2 P aF (z) Mx ¥ ; /x 2 ). (2) 

11 ^ a'=q,g P=G,S,GS 

Here, z G [(3, 1] is the fraction of the Pomeron's light-cone momentum. The first term on the 
right-hand-side of d2J) generates the upper part of the ladder in Fig. [TJa), that is, above /i 2 , 
via the usual (homogeneous) DGLAP evolution with parton-to-parton splitting functions, P aa '- 



3 To be precise, events with a large rapidity gap are described not only by Pomeron exchange, but also 
exchange of secondary Reggeons. From this point of view, LRG events are not strictly diffractive. However, the 
parton densities which originate from secondary Reggeons satisfy the same inhomogeneous evolution equations 
(J2j) . Moreover, even in the inhomogeneous term of |(5J) the Pomeron flux /p(xp;/i 2 ), written in ((3J in terms of 
the inclusive parton densities, also contains some admixture from secondary Reggeons and not just from the 
rightmost (Pomeron) singularity in the complex angular momentum plane. So everywhere the word "diffractive" 
and the index P should be regarded as a symbolic, rather than a literal, notation. 
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The second (inhomogeneous) term describes the cell where the first t-channel single parton is 
generated directly by the Pomeron-to-parton splitting functions, P aP , at a scale /i 2 . In other 
words, these functions describe the transition from the two t-channel partons (that is, the 
Pomeron) to a single t-channel parton. The LO expressions for P aP are given in Ref. [6]. The 
lower parton ladder in Fig. []Ia) plays the role of a perturbative Pomeron flux factor, /p(xp; Li 2 ). 
These flux factors may be calculated from the familiar global parton densities [6]. For example, 
if the two t-channel partons are gluons, then 

(3) 

where the diffractive slope parameter Bp comes from the t integration, and R g from the 
skewedness of the gluon distribution [8] . The notation P = G, S, GS denotes whether these 
two t-channel partons are gluons (P = G) or sea-quarks (P = S), while the interference term is 
denoted by P = GS. We neglect the case where the two t-channel partons are valence quarks, 
assuming that x P <C 1. We will take B D = 6 GeV~ 2 [1], apart from for heavy quark production 
via a direct Pomeron mechanism, Fig. [H(c), where we take Bp = 4 GeV -2 [9,10]. 

Now we come to an important observation. If we neglect the /i 2 dependence of g(x P ,fi 2 ) 
in 0, then / P=G (a;p; /i 2 ) ~ l//i 2 and the inhomogeneous term would be a higher- twist power- 
like correction, which is usually neglected and so collinear factorisation would hold in DDIS. 
However, as x p —>■ 0, the gluon density grows as (/i 2 ) 7 , with anomalous dimension 7 = 0.5, 
while in the HERA domain 7 ~ 0.3. This compensates for the \j [i 2 power-like suppression in 
We therefore cannot neglect the inhomogeneous term in §2§. 

Using only the perturbative framework, that is, just the diagrams in Figs. [T]^a,c), down 
to a very low scale would, of course, be unreliable. However, it is reasonable to assume that 
all the contributions originating from these low scales, /1 2 < /i 2 ,, are parameterised in terms 
of an unknown non-perturbative input, which is schematically presented in Fig. [TJ(b). Unlike 
the unknown input required to analyse inclusive DIS which depends on a single momentum 
fraction, here we have an unknown input which depends on two momentum fractions, x P and 
z. In this non-perturbative domain it is reasonable to assume a Regge behaviour, in which the 
x P dependence is given by the sum of Regge poles (including possible interference terms), and 
that each pole generates input diffractive parton distributions with a freely-parameterised z 
behaviour, where the parameters are to be determined by fitting to the DDIS data. Note that 
for very low x p it would be sufficient to keep just Pomeron exchange. 

Thus in the pQCD approach we consider the DPDFs to be a sum of Pomeron and secondary 
Reggeon contributions, 

a D (xp, z, Q 2 ) = a D < p (xp, z, Q 2 ) + a D - R (x P , z, Q 2 ), (4) 

where we have neglected the possible interference term, as is the usual practice [2] . We param- 
eterise the input quark-singlet and gluon DPDFs for the leading Pomeron pole at a starting 
scale /ip in the form: 

zTP^{x ¥ ,z^l) = /p(x P ) A q z B «(l -z)°\ (5) 
zg D ' v (x P ,z^ 2 ) = f r (xp) A 9 z B s(\-zf^ (6) 



fF=G(x P ; fx 2 ) 
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where the Pomeron flux factor is taken from Regge phenomenology [2], 

^min 

Mxf)=A p / dt e Bpt xl~ 2ar{t \ (7) 

with ofp(t) = ap(0) + a' w t. The secondary Reggeon contributions may be included in a similar 
way. However, in practice, the DDIS data at large Xf are insufficient to constrain the input 
secondary Reggeon parameters. Instead, following HI [2], we fix the z and Q 2 dependence 
of a D ' E (xp, z, Q 2 ) to be the same as the pionic parton distributions, that is, we neglect the 
inhomogeneous terms in the evolution, and fit only an overall normalisation factor cp. This is 
reasonable since the secondary Reggeons are important only for larger xp, where the anomalous 
dimension is not large enough to cancel the power suppression in the corresponding perturbative 
Pomeron flux factors. Moreover, in the perturbative domain, the secondary Reggeons are 
mainly non-singlet quarks which, even at very small xp, have a small (or negative) anomalous 
dimension. This completes the summary of the pQCD procedure which is necessary to describe 
DDIS. 

Now we introduce the Regge factorisation approach and relate it to the above pQCD pro- 
cedure. In the "Regge" description the last term in ([1]) and the inhomogeneous term in ([2]) are 
neglected. The factorising Regge input, ([5]) and (Q, is used in the remaining, now homoge- 
neous, DGLAP evolution equations, ([2]). This assumption is justified for sufficiently large fi^ 
and fixed xp. The smaller the value of xp, the larger the value of /!q has to be to justify the 
pure Regge approach^. What happens if the Regge approach is applied for smaller values of 
/ip? The non- negligible last term in (T5]) is attempting to be mimicked by enlarging the Regge 
input at /2q. As a consequence we obtain the wrong /i 2 dependence of the DPDFs. We can 
readily glimpse what the result will be. Experimentally we measure F^ 3 \ which is dominantly 
determined by the diffractive quark density. Thus, to mimic the /i 2 dependence of the last term 
in ([2]) in a pure Regge fit, we would expect to obtain a larger gluon at relatively large z so as 
to provide a faster growth of the quark density with increasing /i 2 . 

How different are the diffractive parton distributions obtained from fits to the DDIS data 
using the "pQCD" and the simpler "Regge" approaches? To see this, we perform fits to the 
latest HI LRG data [2] for the diffractive reduced cross section . We choose this particular 
data set because HI have already performed a "Regge" parton analysis to precisely these data 
with M x > 2 GeV and Q 2 > = 8.5 GeV 2 . We will use the same kinematic cuts, in order 
to allow a direct comparison with the HI 2006 analysis [2]. We take ap(0) as a free parameter 
in © and fix B F = 5.5 GeV~ 2 and a' ¥ = 0.06 GeV" 2 [1], as in the HI fit. Moreover, the 
secondary Reggeon contributions are also included in a similar way to the HI 2006 analysis [2], 
with a R (0) = 0.5, a' R = 0.3 GeV" 2 and B u = 1.6 GeV" 2 [1], but using the GRV NLO pionic 
parton distributions [11] rather than the Owens LO set [12]. The number of parameters in the 
pQCD and Regge descriptions are the same, since the inhomogeneous term and direct Pomeron 
contributions are calculated explicitly [6], at leading order, and have no free parameters. We 

4 By "pure Regge approach" we mean that we can neglect the inhomogeneous term in ([2]) and describe the 
DPDFs aP by the linear evolution equations starting from input generated by some set of Regge poles. Of 
course, at a large scale it may be necessary to have a more complicated structure in the complex angular 
momentum plane than simply the Pomeron and secondary Reggeon poles. 
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Table 1: Input parameters (corresponding to My = m p ) determined in the Regge and pQCD fits 
to HI LRG data [2] with Q^in = 8-5 GeV 2 . Notice that the A% 2 = 1 errors on the parameters 
are generally smaller in the pQCD fit. 

choose to take the input parameterisations at /i 2 , = 2 GeV 2 , again to keep close to the value 
used in the HI analysis [2], and because the inclusive parton densities which appear in the 
pQCD approach are not well known below 2 GeV 2 where there are no data included in the 
global fits. We use the MRST2004F3 NLO PDFs [13] in the pQCD fits. The inhomogeneous 
evolution is implemented via a modified version of the qcdnum [14] NLO DGLAP program. 

For the Pomeron contribution, we include heavy quarks in the fixed flavour number scheme 
(FFNS), where there are no heavy quark DPDFs, with m c = 1.43 GeV and nit, = 4.3 GeV, 
and use a three-flavour as with Aq^ 3 ' 1 = 407 MeV [13]. For the resolved Pomeron case, 
heavy quarks are generated via photon-gluon fusion at NLO [14] with the renormalisation and 
factorisation scales taken to be (Q 2 + 4ml) where h = c,b. For the secondary Reggeon case, 
the contribution to the diffractive structure functions from heavy quarks is included using the 
zero-mass variable flavour number scheme, and we take the same values of as, m c and rrib as 
in the determination of the pion PDFs [11]. 

The direct Pomeron terms C2,p = CV,p + Cl,p, that is, the second terms on the right-hand- 
side of (JTJ) , are calculated at LO using: 



Ct,f=g = 


(A.27) 


- (A.28), 


Cl,i 


»=g = (A.36), 


(8) 


Ct,¥=S = 


(A.68) 


- (A.69), 


Cli 


»=5 = (A.74), 


(9) 


Ct,v=gs — 


(A.86) 


- (A.87), 


Cl,i 


» =G5 = (A.91). 


(10) 



where the equation numbers refer to the Appendix of Ref. [6]. For light quarks, the lower limit 
of the /x 2 integration in these equations is replaced by /x 2 ,, since the contribution from fi 2 < fi^ 
is included in the input DPDFs taken at a scale /i 2 ,. For CV,p the light-quark contributions 
in the limit fi 2 Q 2 are subtracted since they are already included in the first term of ([T]). 
This subtraction defines a choice of factorisation scheme. There is no such subtraction for C^.p, 
which are purely higher- twist, or for the heavy quark contributions since we are working in the 
FFNS. 

Examples of the diffractive parton densities obtained from these fits to the HI LRG data are 
shown in Fig. [2j During the fits, statistical and systematic experimental errors were added in 
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Q 2 = 25 GeV 2 



Q 2 = 90 GeV 2 




Figure 2: The "Regge" and "pQCD" DPDFs obtained from fitting to the HI LRG data [2] 
with Q^ in = 8.5 GeV 2 compared to the HI 2006 Fit A [2]. All the DPDFs plotted here are 
normalised to My < 1.6 GeV by multiplying by a factor 1.23 relative to My = m p . Only 
the Pomeron contribution is shown and not the secondary Reggeon contributions which are 
negligible at the value of xp = 0.003 chosen here. Also shown are the DPDFs obtained by 
HI [15], using the Regge factorisation approach, from a preliminary combined analysis of their 
inclusive DDIS data [2] with preliminary data on diffractive dijet production [15]. 
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quadrature. An overall factor of 1.23 [1] was applied to the theoretical prediction to account for 
the effect of proton dissociation up to masses My < 1.6 GeV present in the HI LRG data [2]. 
The parameters are given in Table [Hand correspond to My = m p . As is to be expected, there 
are only very small differences between our "Regge" fit and the HI 2006 Fit A0 The value of 
the ap(0) parameter in the pQCD fit is 1.106 ± 0.005. It is a bit higher than the value of the 
intercept 1.08 [5] coming from soft hadron phenomenology, since the input distributions are 
taken at /i 2 , = 2 GeV 2 , and already account for some pQCD contribution which is known to 
be present at this scale. For example, in going from p to J/ip photoproduction, the effective 
Pomeron intercept increases from about 1.09 [16, 17] to 1.20 [9, 10]. Note that a detailed 
comparison of the intercepts would require a study of the correlation between the fitted «p(0) 
and the assumed a' v , whereas here we have fixed a' v = 0.06 GeV~ 2 to agree with the HI 
analysis [2]. 

From Fig. [2] we see that the diffractive quark singlet distribution, S D ' P , is relatively stable 
in going from the pQCD approach to the simplified Regge approach. However, as anticipated, 
the diffractive gluon distribution, g D,p , is greatly suppressed for z > 0.5 in the pQCD approach 
relative to the gluon obtained using the Regge approach. 

The results obtained here are stable to the inclusion of DDIS data from the HI FPS [1], the 
ZEUS leading proton spectrometer (LPS) [18], and the ZEUS diffractive charm reduced cross 
section [19], as well as being stable to including HI LRG data [2] down to Q 2 = 6.5 GeV 2 . At 
an input scale of /i 2 , = 2 GeV 2 , the results are also stable to an extra parameter in the input 
distributions in the form of an additional factor (1 + D a ^fz) and to the use of an alternative 
input parameterisation in terms of Chebyshev polynomials, although this is no longer the case 
if a higher input scale is used such as /i 2 , = 6 GeV 2 . 

When the forthcoming ZEUS LPS and LRG data [20] are released we will perform an analysis 
of the combined data sets, including also the measurements of diffractive dijet production. The 
pQCD diffractive parton densities, a D,p and a D ' M , of the present analysis are publically available^ 
in the form of three-dimensional grids in (xp, z, Q 2 ), together with the corresponding predictions 
for the diffractive structure functions F^ 3 ^, F^ 3 \ F^^' cc an d p^^ ,cc _ The resolved Pomeron, 
direct Pomeron and secondary Reggeon contributions to the diffractive structure functions are 
provided separately in each of these four cases. Measurements of the diffractive charm reduced 
cross section and the fractional contribution of charm to the total diffractive cross section have 
recently been shown to be in excellent agreement with these predictions; see Figs. 12 and 13 of 
Ref. [21]. 

The HI Collaboration have recently presented a preliminary combined fit [15], within the 
Regge factorisation framework, of their inclusive DDIS data [2] with preliminary data on diffrac- 
tive dijet production [15]. The DPDFs from this combined fit, shown in Fig. [2J have a much 
smaller diffractive gluon density at moderate to high z compared to the HI 2006 Fit A. More- 
over, the \ 2 f° r the 190 inclusive DDIS points increases from 158 (HI 2006 Fit A) to 169 (HI 
combined fit) on inclusion of the dijet data. Therefore, the gluon density determined indirectly 



5 These are due to minor details in the two fits such as, for example, the treatment of experimental errors, 
the precise form of the input parameterisation, the choices for as, rn c , mj, the heavy quark renormalisation 
and factorisation scales, and the pion PDFs. However, as seen in Fig. [2j these minor details are not important. 

6 http : //durpdg . dur . ac . uk/hepdata/mrw . html 
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from the inclusive DDIS data, under the assumption of pure DGLAP evolution, is different from 
the gluon density preferred by the dijet data. On the other hand, our pQCD fit naturally has a 
smaller gluon density at moderate to high z than the HI 2006 Fit A, due to the presence of the 
extra inhomogeneous term in the evolution equations. Therefore, the apparent tension between 
the inclusive DDIS and dijet data in the Regge factorisation approach is partly alleviated by 
the inclusion of the perturbative Pomeron terms. As a consequence, we expect that our pQCD 
fit will give a better description of the diffract ive dijet data than the HI 2006 Fit A. 
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